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This document relates to collaborative research between the Complex Systems Research 
Center at the University of New Hampshire and the Ecosystems Center of the Marine Biological 
Laboratory (Woods Hole, MA) under NASA Grant NAGW-714. The authors have organized the 
report as a brief synopsis, graphics, and a collection of published manuscripts which detail key 
methodologies and results. 

The main objective of our research was to characterize biogeochemical cycles at continental 
and global scales in both terrestrial and aquatic ecosystems. This characterization applied to both 
natural ecosystems and those disturbed by human activity. The primary elements of interest were 
carbon and nitrogen and the analysis sought to quantify standing stocks and dynamic cycling 
processes. The translocation of major nutrients from the tenestrial landscape to the atmosphere 
(via trace gases) and to fluvial systems (via leaching, erosional losses, and point source pollution) 
were of particular importance to this study. Our aim was to develop the first generation of Earth 
System Models. 

Our research was organized around the construction and testing of component biogeochemical 
models which treated terrestrial ecosystem processes, aquatic nutrient transport through drainage 
basins and trace gas exchanges at the continental and global scale. A suite of three complementary 
models were defined within this construct. The models were organized to operate at a 1/2 degree 
latitude by longitude level of spatial resolution and to execute at a monthly time step. This 
discretization afforded us the opportunity to understand the dynamics of the biosphere down to 
subregional scales, while simultaneously placing these dynamics into a global context. 

The first model in the suite to be developed was the global drainage basin model, consisting of 
a coupled Water Balance/Water Transport Model (WBM/WTM). This model calculates the water 
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balance by transforming complex patterns of climate into soil moisture, evapotranspiration and 
runoff on single 1/2 degree cells. In addition, it simulates fluvial transport through a networking 
algorithm linking whole systems of grid cells. The next simulation is the Terrestrial Ecosystem 
Model (TEM), a geo-referenced nutrient cycling model for C and N standing stocks and nutrient 
fluxes. The TEM outputs standing stocks of C and N over the year, net primary production and 
respiration. 

Numerous data sets were assembled and analyzed as part of this work. Data required for the 
WBM/WTM and the TEM include climatic drivers for irradiance, temperature, and precipitation, 
and spatial information regarding soil texture, vegetation coverage and elevation. The TEM 
requires a collection of ecosystem-specific datasets for initial calibration and a regional sensitivity 
analysis to remove extraneous input parameters or to restructure the simulation. A total of 42 such 
datasets were compiled and the relevant information catalogued within a computerized database. A 
subset is shown in Table 1. Each data set was processed to conform to our model requirements. 
Point data and spatial data at coarse scales for climate-driving variables were interpolated to the 1/2 
by 1/2 degree latitude by longitude scale using a spherical interpolator developed by Willmott et al. 
(1985), at the University of Delaware. We developed a new map of potential vegetation classes at 
the 1/2 degree scale for the globe (Kicklighter et al., 1994) using Matthews (1983) and extant maps 
at various scales where available. These were compared to potential distributions derived from the 
bioclimatic driving variables from our climate data sets. 

Both the WBM/WTM and the TEM have been successfully implemented in major dryland 
ecosystems at both the continental and global scales. Our work has demonstrated that a relatively 
simple hydrologic model, in conjunction with a suite of datasets at 1/2 degree spatial resolution, 
can generate reasonable hydrodynamic information over both space and time for use in global 
modeling studies. For example, WBM/WTM model results for South America mimic observed 
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datasets for both the pattern of runoff generation and the resultant discharge hydrographs in the 
Amazon river (Vorosmarty et al 1989). The model was subsequently modified to account for 
submonthly soil water and evapotranspiration dynamics (Rastetter et al., 1991) and applied to the 
Zambezi drainage basin in southeastern Africa (Vorosmarty and Moore 1991, Vorosmarty et al 
1991). This work is currently being expanded to cover the major river systems of the globe. 

The TEM has been applied to South America and has produced the most spatially detailed 
estimates of monthly and annual fluxes of CO 2 and net primary production currently available 
(Raich et al., 1991). Application of the TEM to North America revealed that under elevated 
temperature net primary production increases primarily due to the enhanced availability of 
mineralized nitrogen (McGuire et al., 1992, 1993). Globally, the TEM has been used with climate 
data generated by four general circulation models (GCMs) to estimate changes in patterns of net 
primary production over the globe due to a doubling of atmospheric carbon dioxide (Melillo et al., 
1993). 

Early efforts at linking terrestrial and aquatic process models utilized simple accounting 
models both for landscape/mobilization and aquatic transport/retention of nutrients (Gildea et al., 
1986, Vorosmarty et al., 1986) in large rivers (i.e., the Mississippi River). Although this 
approach provides an important tool for quantifying the translocation of materials across the 
continental land mass, it cannot predict the impacts of land use or greenhouse warming. A process 
modeling approach is needed. To this end, we have sought to link the WBM/WTM with the TEM. 
The result will be a mechanistic view of how carbon and nitrogen are first mobilized from the 
terrestrial landscape and transported through fluvial ecosystems with eventual delivery to the 
coastal ocean. We developed a hydrological database of discharge for the 25 worlds largest 
catchments in order to generate global estimates of water and nutrient transport The simulated 
topology for these rivers is shown in Figure 11. 
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With the emergence of Earth Systems Science has come the realization that there exist key 
linkages among the atmosphere, terrestrial biosphere and the worlds oceans. With this realization 
have also come calls for the development of geographically-specific models and associated data 
sets capable of describing the dynamics of ecosystems at continental and global scales. The 
collected works that have emerged from this NASA-funded initiative provide ample evidence that it 
is possible not only to quantify the existing state of key biogeochemical cycles over broad 
geographic domains, but to gain insight into the issue of anthropogenic change. The project has 
produced a first set of component Earth System Models, with a particular emphasis on terrestrial 
and aquatic dynamics. The process-based tools so developed have provided important guidance in 
our current spectrum of research activities, in particular those related to our NASA Earth 
Observing System Interdisciplinary Investigation. 
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FIGURES 


Figure 1. The water balance model showing pools and water transfers. Through a retention 
function for soil water, the model transforms precipitation and potential evapotranspiration into soil 
moisture, evapotranspiration and runoff. Closed triangles represent inputs; open triangles are flux 
determinations made by the model. (From: Vorosmarty et al., 1989). 


Figure 2. Seasonal pattern of runoff generation predicted by WBM/WTM for the 
Amazon/Tocantins basin. A northward movement of relatively high runoff is apparent from 
January through July, followed by low runoff in October throughout the majority of the basin. 
The extreme northwestern portion of the river system contributes high runoff throughout the year. 


Figure 3. Observed and predicted mean annual runoff for South America. The observed 
dataset is from Korzoun et al (1977). The predictions were made using the global WBM 
developed under this research. (From: Vorosmarty et al. 1989) 


Figure 4. The terrestrial ecosystem model (TEM). The state variables are: carbon in 
vegetation (Cy); nitrogen in vegetation (Ny); organic carbon in soils and detritus (C<$); organic 
nitrogen in soils and detritus (Ng); and available soil inorganic N (N^y). Arrows show carbon 
and nitrogen fluxes: GPP, gross primary productivity; R^, autotrophic respiration; Rjj, 
heterotrophic respiration; L c , litterfall C; L N , litterfall N; NUPTAKE, N uptake by vegetation; 
NETNMIN, net N mineralization of soil organic N; NINPUT, N inputs from outside the 
ecosystem; and NLOST, N losses from the ecosystem. (From: McGuire et al. 1992) 


Figure 5. Potential net primary productivity (NPP) in South America (as carbon) for the 
months of January, April, July and October. Values are g-m'^-mo'* of carbon. Negative NPP 
values indicate autotrophic respiration exceeded gross primary production during that month. The 
blockiness in these figures is due to the poor spatial resolution of the cloudiness data set used to 
estimate the monthly irradiance of photosynthetically active radiation. (From: Raich etal. 1991) 
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Figure 6. (top) Annual net primary productivity of potential vegetation in North America as 
determined by the terrestrial ecosystem model, (bottom). Annual net nitrogen mineralization of 
potential vegetation in North America as determined by the terrestrial ecosystem model. (From: 
McGuire et al. 1992) 

Figure 7. Response of net primary productivity (NPP) and net nitrogen mineralization 
(NETNMEN) for temperate mixed forests of North America in an elevated-temperature experiment 
Monthly temperature was increased 2° for the terrestrial ecosystem model (TEM) or the water 
balance model (WBM), which provides climatic inputs to TEM. Blank and diagonally-lined bars 
correspond to the NPP response when TEM was run with the carbon cycle uncoupled or coupled 
to the nitrogen cycle, respectively. The crosshatched bar is the NETNMIN response with the 
carbon cycle coupled to the nitrogen cycle. (From: McGuire et al. 1993) 


Figure 8. Potential natural vegetation used by the TEM to estimate annual NPP for the global 
terrestrial biosphere. (From: Melillo et al. 1993) 


Figure 9. Annual NPP estimated by the TEM for the global terrestrial biosphere. (From: 
Melillo et al. 1993) 


Figure 10. Percent difference in annual NPP between contemporary climate at 312.5 
p.p.m.v. CC >2 and the various GCM climates at 625.0 p.p.m.v. CO 2 , as predicted by the TEM for 
eastern North America (left) and southeast Asia (right). Estimates were not made for open water or 
wetland ecosystems. (From: Melillo et al. 1993) 


Figure 11. Network topology for the 25 largest rivers of the globe. 
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TABLE I Sample data holdings at UNH. 


TEMPORAL 


DATATYPE 

SOURCE 

SPATIAL SCALE 

RESOLUTION 

Physical Parameters 




Coastlines and Political Boundaries 

World Data Bank D (CIA/NTIS) 

Vector, rasterized to 0.5° x 0.5° 

N/A 

Elevation and Terrain Characteristics 

NAVY (1984) 

10 x 10 

N/A 


ET0P05 (1988) 

5' x 5' 


Development Indicators 




Methane Emission from Animals 

Matthews and Fung (1987) 

1.0° x 1.0° 

N/A 

Population 

World Data Bank 

By country 


Urban Development 

NAVY (1984) 

10 x 10 

N/A 

Climate 




Air Temperature 

Legates and Willmott (1988) 

0.5° x 0.5° 

monthly 


Spangler and Jenne (1988) 

Available stations 

monthly to 1988 

Precipitation 

Legates and Willmott (1988) 

0.5° x 0.5° 

monthly 


Spangler and Jenne, (1988) 

Available stations 

monthly 


Nicholson (1988) 

Available stations (Africa only) 

monthly 


NOAA 

Available stations (Africa only) 

monthly, 10-day 

Cloudiness, Net Irradiation, 

Hahn et al. (1988) 

varies, 2.5°-30.0° 

quarterly 

and PAR 

3DNEPH/US AFETAC 

25 nautical miles 

daily, 1970-77 


RTNEPH/USAFETAC 

25 nautical miles 

daily, from 1973 

Radiation, Solar and Thermal 

NCAR (various GCM results) 

Grid format, various resolutions 

monthly 

Wind, u- and v-components 

GRID/FAO (1986) 

Vector format 

annual 


NCAR (GCM results, 1989) 

Grid format, various resolutions 

monthly 


HIRAS/USAFETAC 

2.5° x 2.5° 

daily 

Soils 

F AO/C SRC (1985) 

0.5° x 0.5° 

N/A 


USSR Inst of Geogr./CSRC 

Vector, rasterized to 0,5° x 0.5° 

N/A 

Vegetation 




Potential vegetation and Cultivation 

Matthews (1983) 

1.0° x 1.0° 

N/A 

Land Cover and Land Use 

White (1983) 

Vector, rasterized to 0.5° x 0.5° 

N/A 

Holdridge Life Zones 

Emanuel et al. (1985) 

0.5° x 0.5° 

N/A 

Contemporary Land Use 

USSR Inst. of Geogr./CSRC 

Vector, rasterized to 0.5° x 0.5° 

m 

Potential Ecozone Classification 

MBL/CSRC (1990) 

0.5° x 0.5° 

IN /A 

Ecosystem Complexes 

Olson et al. (1985) 

0.5° x 03° 

N/A 

Hydrological Data 




Discharge hydrographs 

UNESCO (various years) 

Vector format 

monthly 

Direction of flow 

Operational Nav. Charts 

Vector format 

N/A 

Runoff 

WMO/Global Runoff Center 

Available stations 

monthly 

Watersheds 

GRID/FAO (1986) 

Vector format 

N/A 
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